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ABSTRACT 

We analyze photometry and spectroscopy of a sample of 63 clusters at 0.3 < z < 0.9 drawn from 
the Las Campanas Distant Cluster Survey to empirically constrain models of cluster galaxy evolution. 
Specifically, 1) by combining /—band photometry of 44 of our clusters with that of 19 clusters from 
the literature (Aragon-Salamanca et al. ; Smail et al. ; Stanford et al. ) we parameterize the redshift 
dependence of M*[ in the observed frame as M/ = (-21.74±0.12) - (0.88±0.24)z-5 log/i (rms deviation 
= 0.34) for 0.3 < z < 0.9 (O = 0.2, n A = 0)); 2) by combining 30 of our clusters and 14 clusters 
from the literature (Aragon-Salamanca et al. ; Smail et al. ; Stanford et al. ) with V and / data we 
parameterize the redshift dependence of the V — I color of the E/SO red sequence in the observed frames 
as V-I = (-0.24±0.28) + (7.42±1.03)z-(4.61±0.91)z 2 (rms deviation = 0.16) for 0.3 < z < 0.9; and 3) 
by combining 13 of our clusters with 15 clusters from the literature (Aragon-Salamanca et al. ; Stanford 
et al. ) with / and K 1 data we parameterize the redshift dependence of the I—K' color of the E/SO red 
sequence in the observed frames as I-K' = (0.66±0.65) + (9.50±3.72)z-(14.72±7.01)z 2 +(8.72±4.29)z 3 
(rms deviation = 0.18) for 0.3 < z < 0.9. Using the peak surface brightness of the cluster detection, E, 
as a proxy for cluster mass (Gonzalez et al. ), we find no correlation between S and MJ or the location 
of the red envelope in V — I. We suggest that these observations can be explained with a model in 
which luminous early type galaxies (or more precisely the progenitors of current day luminous early 
type galaxies) form the bulk of their stellar populations at high redshifts (.> 5) and in which many of 
these galaxies, if not all, accrete mass either in the form of evolved stellar populations or gas that causes 
only a short term episode of star formation at lower redshifts (1.5 < z < 2). Our data are too crude to 
reach conclusions regarding the evolutionary state of any particular cluster or to investigate whether the 
morphological evolution of galaxies matches the simple scenario we discuss, but the statistical nature of 
this study suggests that the observed evolutionary trends are universal in massive clusters. 

Subject headings: galaxies: clusters: general — galaxies: evolution — galaxies: formation — galaxies: 
photometry — surveys 



1. INTRODUCTION 

High redshift clusters provide an excellent laboratory 
with which to study the evolution of galaxies in dense envi- 
ronments. By investing a modest amount of telescope time 
obtaining spectra of a handful of galaxies, thereby verify- 
ing the reality of the cluster and ascertaining its redshift, 
one is rewarded with the knowledge of distances to tens or 
even hundreds of galaxies. The spectral and photometric 
properties of these galaxies can then be used to constrain 
evolutionary models, whose differences usually become ap- 
preciable at high redshift, z > 0.5 (cf. Bruzual & Chariot 
1993; Kauffmann et al. 1998). Until recently, such efforts 
have been hampered by the small number of known high 
redshift clusters. Published work in this field relies upon 
the examination of the same ~dozen clusters, drawn pri- 
marily from two optical surveys. Both of these surveys 
(Gunn, Hoessel, & Oke 1986, hereafter GHO; Couch et al. 
1991, hereafter C91) selected objects on the basis of an 
overdensity of resolved objects. GHO's efforts yielded 8 



clusters with z > 0.55, while C91 yielded 3. 

Because the number of known high redshift clusters is 
small, they have been studied in great detail by many au- 
thors. In particular, with the advent of the Keck telescopes 
and the Hubble Space Telescope (HST), detailed informa- 
tion regarding the cluster galaxies' star formation rates 
and morphologies have complemented the more traditional 
ultraviolet, optical, and infrared photometry. As a result, 
a comprehensive picture of the formation and evolution 
of cluster galaxies has begun to emerge in which many, if 
not all, galaxies have experienced some recent star forma- 
tion. The most active class of galaxies was first noted by 
Butcher & Oemler (1984, B-O). They studied 33 clusters 
of galaxies with 0.003< z <0.54 and discovered that the 
fraction of blue galaxies increased markedly with redshift. 
They interpret the rise in blue fraction as evidence for in- 
creased star formation in the past and this interpretation 
has subsequently been confirmed by spectroscopic stud- 
ies (Dressier & Gunn 1982, 1983, 1992; Couch & Sharpies 
1987; Couch et al. 1994; Abraham et al. 1996; Fisher et 
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al. 1998; Poggianti et al. 1999). Furthermore, morpho- 
logical classifications using high resolution HST imaging 
show that the bulk of these star-bursting galaxies are disk 
dominated systems, a significant fraction of which show 
signs of interactions or mergers (Couch et al. 1994, 1998; 
Dressier et al. 1994, 1999; Oemlcr, Dressier, & Butcher 
1997; Poggianti et al. 1999). 

Cohabitating with these active galaxies, however, is a 
seemingly quiescent population. These galaxies are the 
cluster's brightest and reddest, which lie along a very nar- 
row locus in the cluster color-magnitude (CM) diagram re- 
ferred to as the red envelope. Detailed morphological and 
spectroscopic studies of clusters at high redshift confirm 
that these galaxies are a subset of the cluster's elliptical 
population (e.g. de Propris et al. 1999; van Dokkum et 
al. 1999). Evolution of the zeropoint, slope, and scatter 
about cluster CM relations all suggest that these ellip- 
ticals formed at a very high redshift (z J> 5) and have 
experienced little or no recent star formation (see Aragon- 
Salamanca et al. 1993, hereafter A93; Rakos & Schombert 
1995, hereafter RS; Ellis et al. 1997; Stanford et al. 1998, 
hereafter SED98; Jones et al. 2000), while work on the 
evolution of the fundamental plane (Kelson et al. 2000) 
and on the evolution of absorption line indicics (Bender, 
Ziegler, & Bruzual 1996; Kelson et al. 2001) support this 
conclusion by placing a lower redshift limit of ~ 2.5 for the 
formation of the bulk of the stellar populations. Studies 
of the optical and infrared luminosity functions of cluster 
galaxies suggest a more moderate formation redshift for 
cluster ellipticals, z ;> 2 (Smail et al. 1997, hereafter S97; 
de Propris et al. 1999). Despite this general consensus, the 
non-uniqueness of spectral synthesis models, against which 
the majority of these observations are compared, allows for 
the possibility that at least some of these galaxies are not 
completely quiescent. In fact, recent work by van Dokkum 
et al. (1999) suggests that a large fraction of these galax- 
ies in the higher redshift (~ 1) clusters are not relaxed 
systems. These authors studied MS1054-03 at z = 0.83 
using spectroscopy and HST high resolution imaging and 
found that 17% of the L>L* cluster population is experi- 
encing a major merger, most of which will probably evolve 
into luminous (2L») elliptical galaxies. They estimate that 
~50% of the present-day cluster ellipticals were assembled 
in mergers since z ~ 1, assuming that the galaxy popu- 
lation of MS1054-03 is representative. While most of the 
mergers involve red early-type galaxies with no detected 
[Oil] A3727 emission, some do exhibit enhanced Balmer 
absorption indicative of a modest recent starburst. Con- 
sequently, the question of the evolutionary history of the 
brightest cluster ellipticals remains open. 

As mentioned previously, all of the above work is ham- 
pered by the lack of statistically significant samples. Gen- 
eral conclusions based on small samples rely critically on 
the assumption that this handful of clusters is represen- 
tative of the entire high redshift cluster population. This 
assumption is particularly worrisome when (1) it is viewed 
in the context of a hierarchical clustering scenario which 
holds that structures of different mass scales form at dif- 
ferent times and consequently have experienced different 
evolutionary histories, and (2) one realizes that the cluster 
properties play a significant role in their original identifi- 
cation. The need for a large, well-defined, high redshift 



cluster sample is evident. 

In an attempt to resolve issues raised by small sample 
size, Postman et al. (1996) conducted the Palomar Dis- 
tant Cluster Survey. This survey consists of optical/near 
IR imaging of a ^5 deg 2 area using the 4-shooter CCD 
camera on the Palomar 5m telescope. Their cluster find- 
ing technique was more sophisticated than previous opti- 
cal surveys by incorporating magnitude and color filters 
to minimize the contamination from line-of-sight projec- 
tions. Their method relies upon resolving distant galaxies 
in two colors and thus requires the use of a large telescope 
under photometric conditions with good seeing, which lim- 
its the areal coverage of the survey. Their efforts yielded 
a large, well-defined cluster catalog that contains 35 clus- 
ters with estimated redshifts >0.5, plus 25 clusters from an 
inhomogeneous sample with estimated redshifts >0.5 (the 
PDCS includes clusters identified previously by GHO). Us- 
ing these clusters to study galaxy evolution, their results 
to date provide more statistically significant confirmation 
of the trends observed with smaller samples. In particular, 
Lubin (1996) examined the location of the red envelope in 
the HST filters V 4 ~h {V 4 = F555W and I 4 = F785LP) and 
found that the data are most consistent with the passive 
evolution of an elliptical galaxy that formed at high red- 
shift. She also measured the fraction of blue galaxies, ff,, 
and found an increase in f;, with redshift, consistent with 
a high redshift extension of the Butcher-Ocmler effect. 

An alternative, efficient method for finding high red- 
shift clusters was developed by Dalcanton (1996). This 
technique capitalizes on the idea that the total flux from 
a distant cluster in shallow images is dominated by the 
flux from unresolved sources. If one has an intrinsically 
uniform image, the flux from the unresolved sources pro- 
duces a detectable low surface brightness fluctuation. Be- 
cause the technique does not rely upon an overdensity of 
resolved objects, relatively shallow (but highly uniform) 
exposures are sufficient to identify distant clusters and a 
much larger area on the sky can be surveyed. Dalcanton 
et al. (1997) verified the ability of this technique to detect 
LSB objects to the required central surface brightness us- 
ing existing scans from a ~ 17.5 deg 2 Palomar 5m survey 
(originally taken as part of the Palomar Transit Grism Sur- 
vey for high redshift quasars; Schneider et al. 1994). This 
survey identified more than 50 northern hemisphere clus- 
ter candidates. In 1995, the authors conducted a 130 deg 2 
drift-scan survey of the southern sky at the Las Campanas 
Observatory. The reductions and object identification is 
complete and the final cluster catalog contains over 1000 
cluster candidates (Gonzalez et al. 2001). Not only does 
this catalog result in a tremendous increase in the number 
of known clusters at these redshifts, but because our clus- 
ter identification criteria is independent of those utilized 
in previous surveys, this catalog provides an independent, 
well-defined sample with which to compare the results of 
more traditional surveys. 

In this paper, we present the first study of galaxy evolu- 
tion using a subset of the Las Campanas Distant Cluster 
Survey. In addition to 54 southern clusters, we include 
9 northern clusters drawn from the Palomar survey. The 
combined sample consists of 63 clusters, 17 of which are 
spectroscopically confirmed and 46 of which have photo- 
metrically estimated redshifts. We obtain deep follow-up 
imaging in V, /, and K' , although not every cluster is 
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imaged in every band. We use these data, along with that 
of 19 clusters from the literature, to classify cluster candi- 
dates, photometrically estimate their redshifts, and study 
the luminosity and color evolution of cluster galaxies. In 
§2 we briefly describe the survey, including cluster can- 
didate identification, and our follow-up observations. In 
§3 we present our three photometric rcdshift estimation 
techniques. In §4 we utilize our clusters with spectro- 
scopic redshifts and those with photometric redshifts to 
trace the evolution of cluster galaxies in luminosity and 
color. Specifically, we examine the evolution of M*j and 
the location of the red envelope in V — I and I—K' . We 
investigate cluster selection biases and their effect on our 
empirically measured galaxy evolution in §4.4. In §4.5 we 
combine our results with those from the recent literature 
and suggest a picture of cluster galaxy formation and evo- 
lution. Finally, we summarize our results and conclusions 
in §6. 

2. THE DATA 

Our data originate from a variety of telescopes and in- 
struments. The candidate galaxy clusters are identified 
using drift scan images and techniques described briefly 
below for context, but in full detail by Gonzalez et al. 2001. 
We then describe the follow-up spectroscopy and photom- 
etry of subsamplcs of candidate galaxy clusters that con- 
stitute the core dataset for this study. 

2.1. Cluster Identification 

Our candidate clusters are drawn from two optical drift 
scan surveys. The first is a ~ 17.5 deg 2 Palomar 5m sur- 
vey using the 4-shooter camera (Gunn et al. 1987) and 
the F555W "wide V" filter described by Dalcanton et al. 
(1997). Dalcanton et al. made use of existing scans (orig- 
inally taken as part of the Palomar Transit Grism Survey 
for high redshift quasars; Schneider et al. 1994) to ver- 
ify the low surface brightness (LSB) technique for detect- 
ing LSB objects to the required central surface brightness 
first proposed in detail by Dalcanton (1996). This survey 
identified more than 50 northern hemisphere cluster candi- 
dates. The second is a 130 deg 2 survey of the southern sky 
conducted at the Las Campanas Observatory lm telescope 
over 10 nights by the authors using the Great Circle Cam- 
era (Zaritsky, Shectman, & Bredthauer 1996) and a wide 
filter, denoted as W, that covers wavelengths from 4500 A 
to 7500 A. This survey was specifically designed to study 
LSB galaxies and high redshift galaxy clusters. Scans are 
staggered in declination by half of the field-of-view (24 ar- 
cmin) so that each section of the survey is observed twice, 
but with a different region of the detector. Exposure times 
are limited by the sidereal rate, which at the declination 
of the survey (—10° to —12°) is ^97 sec. The pixel scale is 
0.7 arcsec/pixel with a typical seeing of about 1.5 arcsec. 
The full region surveyed measures ~ 85° x 1.5° and is cm- 
bedded within the area covered by the LCRS (Shectman 
et al. 1996). The result of iterative flatfielding, masking, 
filtering, and object detection is a catalog of ~ 1000 cluster 
candidates (a full description of the method and catalog 
are provided by Gonzalez et al. 2001). 

The focus of this study is the photometry and spec- 
troscopy of subsamplcs of candidate clusters obtained to 
test cluster identification algorithms and study the clus- 
ters' galaxy populations. Keck spectroscopy, optical (V- 



and /-band) photometry, and IR (if '-band) photometry 
are presented, although not all types of data are available 
for every candidate. Table 1 contains information only 
about those candidates that we determine to be bonafide 
clusters, rather than all observed candidates. The name 
of the cluster is listed in the first column. RA and DEC 
(2000.0) are in the next column, where the center of the 
cluster is taken to coincide with the centroid of the low 
surface brightness fluctuation. The last six columns con- 
tain the exposure time and seeing for each of the three 
photometric bands, V, I, and K' . In summary, 97% of 
our clusters are imaged in I, 54% in V, 33% in K' , and 
29% have spectroscopy. Some clusters were imaged twice 
in various photometric bands. Because objects were ob- 
served on different runs (sometimes with a different cam- 
era and/or telescope) and calibrated using different sets of 
standard stars, they often have two completely indepen- 
dent sets of photometry. We use both sets independently 
in the calibration of our photometric redshift estimators 
and subsequent analysis. 

These data provide a set of confirmed clusters with 
which to develop a set of quantitative statistical descrip- 
tors to use for cluster selection and a statistical sample of 
clusters with which to study galaxy evolution in clusters. 
Because we were testing our cluster selection algorithms, 
the cluster subsample is heterogeneous. It includes both 
some of our best and some of our worst candidates. Be- 
cause our selection criteria evolved during this study, some 
of the confirmed candidates presented here will not be in 
our final cluster catalog. No biases in galaxy properties, 
except those that possibly correlate with cluster proper- 
ties such as concentration and richness (parameters that 
factor into the cluster selection), are expected. This issue 
will be discussed in more detail in §4. 

2.2. Spectroscopy 

Two critical cluster parameters for any subsequent study 
of our catalog are the cluster redshifts and the masses (or 
alternatively to mass, richness, X-ray luminosity, or ve- 
locity dispersion, albeit with different caveats). With a 
catalog of ~ 1000 cluster candidates, it is not feasible to 
obtain spectroscopy or deep follow-up photometry for ev- 
ery candidate. Consequently, it is imperative to develop 
efficient observational estimators of these parameters. We 
present and evaluate several redshift estimators. 

In December 1995, we began an observing program that 
combines long-slit spectroscopy obtained using the Keck 
telescopes and optical and infrared photometry obtained 
using the Las Campanas 2.5m and lm telescopes, the Palo- 
mar 1.5m, and the Lick 3m. Our basic goal is to identify 
and calibrate photometric redshift estimators. These data 
have two additional functions: 1) they enable us to confirm 
candidate clusters and tune our selection criteria, and 2) 
they enable us to study the clusters themselves and their 
galaxy populations. 

Spectra were obtained using the Low-Resolution Imag- 
ing Spectrograph (Oke et al. 1995) on the Keck I telescope 
on 1995 December 20-21, 1997 March 14-15, 1998 April 4- 
5, and 1999 March 20-21. Due to weather and technical 
problems ~ 40% of the time was usable. The deep imag- 
ing necessary to construct multi-object masks was gener- 
ally not available so we used a single long slit with a 600 



4 



Cluster Galaxy Evolution and the LCDCS 



line mm grating. The slit was aligned, using the guider 
image, to include as many individual galaxies as possible 
(generally 2 to 3 with m^ <; 22) and to lie across the posi- 
tion of the LSB feature detected in the scans. As an inter- 
esting aside, the Keck guider images are deeper than the 
drift scan images from which the candidates are selected. 
Typical spectroscopic exposure times were 30 minutes to 
one hour. Preliminary results from this work are discussed 
by Zaritsky et al. (1997). 

The reduction of these spectra is standard. Images are 
rectified and calibrated using calibration-lamp exposures 
and night sky lines (see Kelson et al. 1997 for details). Ve- 
locities are measured using a cross-correlation technique 
or by measuring centroids of emission lines. The rms red- 
shift difference between redshifts measured from different 
emission lines in the same spectra is 0.0033 (typically less 
than 0.0005). In cases where only one emission line is 
observed, we attribute it to H/3. Little ambiguity exists 
between H/3 and [Oil] because the [Oil] emission line is a 
resolved doublet (3726 and 3729 A). 

The distribution of galaxy redshifts is shown in Figure [l] 
for all 28 fields observed. We consider a candidate object 
to be a real cluster if at least three galaxies cluster within 
1000 km s _1 of their mean redshift. A valid concern is 
that "successful" candidates in this scheme may arise ran- 
domly in lines-of-sight through large-scale structures. If 
random large scale structures appeared in 19 of 28 fields 
(60%), then we would expect to find two coherent struc- 
tures in 13 of our fields. Only one such field exists, and so 
it is evident that random redshift associations are not the 
dominant cause of the redshift associations. However, to 
test whether random redshift associations are significant, 
we use published Keck LRIS spectroscopy (Guzman et al. 
1997; Lowentha! et al. , 1997; Phillips et al. 1997; Vogt 
et al. 1997; Cohen et al. 2000) out to z = 1.1 to estimate 
the likelihood of chance redshift associations. The use of 
data from the same instrument as used in this study en- 
sures that the data are of comparable quality and have 
similar systematic problems (sky lines, fringing) over the 
same redshift range. The joint catalog contains coordinate 
positions and redshifts for 686 galaxies at 0.3 < z < 1.1 (a 
redshift range chosen to match that of our sample). We 
then randomly select fields with the same distribution of 
galaxies per field as seen in Figure [l]. The number of pairs 
of galaxies within 1000 km sec^ 1 predicted for our sam- 
ple from 10,000 simulations is 9.9 (we observe 7, which is 
within the la Poisson fluctuation). The number of triplets 
predicted is 1.3 (consistent with the one system in which 
two structures are seen, although this is not necessarily 
the spurious system). Only four of our successful 19 cluster 
candidate have the minimum three members. The number 
of quadruples predicted in this sample is 0.2 (we observe 
15). We conclude that systems with three or more galaxies 
within 1000 km sec -1 are overwhelmingly not the result of 
random lines of sight through large-scale structures. This 
conclusion will later be confirmed by the independent cor- 
respondence of the magnitude of the brightest galaxy in 
each cluster candidate, the mean spectroscopic redshift of 
the galaxies in the overdensity, and the concentration of 
galaxies at the position of the original LSB detection. We 
conclude that our success rate for identifying true clus- 
ters is J> 60%, and that perhaps it is substantially greater. 



This success rate will be confirmed independently by other 
means below. 

2.3. Optical Photometry 

V- and /-band images, which straddle the Ca II H & 
K break for galaxies at redshifts of 0.3 <J z <J 0.9, pro- 
vide a measure of the stellar populations in these galaxies. 
Because the majority of our photometry was completed 
before the cluster finding algorithm was fine-tuned (in- 
deed it was necessary to help train the algorithm) , we im- 
aged some candidates only in /, to provide more images 
of potential clusters at the expense of color information. 
Optical photometry in the Johnson y-band and Cousins 
/-band for the northern hemisphere clusters was obtained 
at the Palomar 1.5 m telescope during 1996 May 9-14. 
Conditions during this run were rather poor due to inter- 
mittent cirrus with seeing ranging from 2" to 3". Typical 
exposure times were between 0.5 and 1.5 hr. The south- 
ern hemisphere clusters were imaged in V and I at the 
Las Campanas Observatory using the 1 m (1997 February 
26-28 & March 1-7, 1998 March 27-31 & April 1) and 2.5 
m (1996 March 24-26, 1997 April 11-15, and 1998 March 
28-29) telescopes. Conditions during these runs varied, 
but were predominantly photometric with seeing typically 
~1" (often sub-arcsecond) on the 2.5 m and ~1.3" on the 
1 m. Exposure times ranged between 0.75 and 1.5 hr on 
the 1 m and 0.3 and 1 hr on the 2.5 m. 

Data reduction is once again fairly standard. Individual 
frames are bias subtracted, flatfielded with either twilight 
or dome flats, and combined with appropriate high and 
low sigma clipping to remove cosmic rays and dead pixels. 
Calibration is done using Landolt's (1983, 1992) standard 
fields. All fields observed at Las Campanas were observed 
at least once during photometric conditions. Although the 
nights during the Palomar run were generally not photo- 
metric and bootstrapping photometric exposures are not 
available, the solutions for the standards taken during the 
run show only modest scatter (0.06 mag) that suggests 
the conditions were reasonable. Nevertheless, those data 
should be viewed with caution. 

Galaxy magnitudes are measured using SExtractor's 
(Bertin & Arnouts 1996) "best" total magnitudes. For 
uncrowded objects, this is the automatic aperture magni- 
tude, which misses less than 6% of the light from a galaxy. 
The corrected isophotal magnitude, which misses less than 
15% of the galaxy's light but is more robust than the au- 
tomatic aperture magnitude, is used for crowded objects. 
Although the internal uncertainties in the SExtractor pho- 
tometry almost always dominate, we include the uncer- 
tainty from the photometric solutions in determining the 
uncertainty in the final magnitudes. Galaxy magnitudes 
are corrected for extinction using the dust IR emission 
maps of Schlegel et al. (1998), but are not K-corrected. 

2.4. IR Photometry 

We image a subset of clusters in the infrared K 1 band 
(2.16/im) to probe the older, underlying stellar popu- 
lations. A related advantage is that K-corrections are 
smaller than for the optical passbands, simplifying trans- 
lations to the rest frame magnitudes. The principal disad- 
vantage is decreased efficiency due to the bright sky back- 
ground at K'. 
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Once again, we used several telescopes to collect the 
data. The northern hemisphere clusters and a few south- 
ern hemisphere clusters were observed at the Palomar 1.5 
m telescope (1996 January 7-12), and the Lick 3 m tele- 
scope (1996 April 26-29 and 1998 March 6-9). A majority 
of the nights were photometric, although the seeing was 
rather poor (1" to 2"). Exposure times, on object, were 
typically 1 to 2 hr at Palomar and 0.6 to 1 hr at Lick. 
We observed the bulk of the southern clusters at the Las 
Campanas 2.5 m telescope during 1996 March 31 and 1996 
April 1-2. Every night of this run was photometric with 
excellent seeing (0.6" to 1"). Typical exposure times were 
0.5 to 0.75 hr on object. The telescope was always dithered 
to nearby blank fields at frequent intervals to monitor the 
rapidly changing sky. 

Data reduction follows the standard IR procedure. We 
subtract the dark current from individual images, flat- 
field using either twilight or dome flats, and subtract the 
sky. Images are aligned and combined with an appro- 
priate sigma clipping to remove cosmic rays, and hot or 
dead pixels that were not masked. Every cluster with IR 
photometry was observed at least once during photomet- 
ric conditions. Calibration was performed using the HST 
faint near-IR standards (Persson et al. 1998). The photo- 
metric solutions for the standards taken during the 1996 
Lick run show large scatter (0.1 mag) cither because the 
nights were not photometric or because only a few obser- 
vations of standard stars were made. Consequently, the 
photometric solutions are not sufficiently well determined 
that one or two outliers can be confidently removed. K' 
magnitudes were measured using SExtractor's (Bertin & 
Arnouts 1996) "best" total magnitudes (see above for a 
description). As with the optical photometry, the uncer- 
tainty in the final magnitudes is determined by folding in 
the error in the photometric solutions with SExtractor's 
internal photometric uncertainty. Again, the galaxy mag- 
nitudes are corrected for extinction (Schlcgcl et al. 1998), 
but are not K-corrected. 

3. PHOTOMETRIC REDSHIFTS 

The first step in any potential study of the galaxy cluster 
candidates is to obtain rcdshift measurements. Obtaining 
spectroscopic redshifts for a sample as large as ours is not 
feasible, so we must construct a reliable, but efficient pho- 
tometric redshift technique. We develop and test three 
techniques. First, we use the observed magnitude of the 
brightest cluster galaxy (BCG) measured on our drift scan 
images to obtain zbcg- Second, we use the /—band lumi- 
nosity function measured from the deeper /-band follow-up 
observations to obtain zi um . Third, we use the location of 
the "red envelope" in the cluster color-magnitude (CM) 
diagram measured using the V- and /-band follow-up ob- 
servations to obtain z co i . We empirically calibrate all three 
methods using a subset of our clusters with spectroscopic 
redshifts. In addition, we include 19 clusters taken from 
the literature (A93; S97; SED98) in the calibration of zi um 
and z co i . An empirical calibration is advantageous because 
the estimated redshifts are independent of assumptions of 
cosmology or galaxy evolution. Recall that several of our 
clusters were observed twice in various bands on differ- 
ent runs. Because they have independent calibration and 

1 We cannot reliably measure the cluster surface brightness of about 



photometry, we use both data sets in our calibration. 

The methods vary in robustness and efficiency. The 
most observationally efficient is zbcg because no follow- 
up observations arc necessary (aside from those required 
to calibrate the method). We can measure the magnitude 
of the BCG directly from the drift scan images out to 
z ~ 0.8. Consequently, this is the only redshift indicator 
that we have for all of our cluster candidates with z <J 0.8. 
However, it would appear to be the least robust option be- 
cause it depends on the proper identification of one galaxy 
and on the homogeneity of BCG luminosities at a given 
redshift. Fitting the galaxy luminosity function of each 
cluster and using as a standard candle circumvents 
issues related to small number statistics and the proper 
identification of the BCG. Unfortunately, other difficulties 
arise. The cluster galaxy catalog must be corrected for 
contamination from interlopers by statistically subtract- 
ing the background galaxies, which introduces a significant 
source of uncertainty. Furthermore, the method is intrinsi- 
cally less efficient because it requires follow-up imaging of 
each cluster. The third method, which utilizes the location 
of the red envelope, does not require background subtrac- 
tion because the red ellipticals in the cluster are typically 
the reddest objects in the field, but requires follow-up ob- 
servations in two separate filters. Although no method ap- 
pears perfect, they provide multiple cross checks against 
each other, and against the spectroscopic redshifts. As we 
show next, they all work remarkably well. 

3.1. Estimating Redshifts Using BCGs 

BCGs are excellent standard candles at low redshift, 
z <j0.05, with absolute magnitude dispersion of 0.30 mag 
in V and 0.24 mag in R, when corrections based on envi- 
ronment are applied (e.g. Sandage & Hardy 1973; Hoes- 
sel, Gunn & Thuan 1980; Hoessel 1980; Schneider, Gunn, 
& Hoessel 1983a, 1983b; Lauer & Postman 1992; Post- 
man & Lauer 1995). More recently, their near constancy 
in luminosity has been shown to extend to z ~1 in the 
if-band (A93; Aragon-Salamanca et al. 1998; Collins & 
Mann 1998). Because of this stability, their large intrinsic 
luminosity, and the fact that we can measure their mag- 
nitudes directly from the drift scan images, they are the 
most accessible photometric redshift indicator available to 
us. 

We define the BCG to be the brightest galaxy within 
15" of the candidate cluster center. A fundamental con- 
cern is that BCGs may not necessarily lie coincident with 
the galaxy luminosity centroid measured by our LSB tech- 
nique. Although this is a valid concern, the overall suc- 
cess of the method argues that while any particular cluster 
may be affected, it is not a general problem. (This issue 
is discussed in more detail by Nelson et al. 2001.) To ob- 
tain the BCG magnitude, we measure the light within a 
5" aperture using SExtractor. A second valid concern is 
that a fixed angular aperture, which will correspond to 
different physical scales at different redshifts, is used to 
determine the magnitude. However, because we calibrate 
the magnitude-redshift relation empirically, we are only 
assuming that at each redshift the population is uniform 
(any bias introduced by using different metric apertures at 
different redshifts will calibrate itself out of the relation). 

1/3 of the cluster candidates identified in our earliest analysis because 
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We use a subset of 17 clusters with spectroscopic red- 
shifts from our survey to calibrate the magnitude-redshift 
relation, shown in the upper panel of Figure ^. We model 
this relationship with a function that includes both lumi- 
nosity evolution and a second-order correction term based 
on the surface brightness (£) of the cluster in the detection 
data (see Gonzalez et al. 2001) to constrain the function to 
have the generally understood behavior of distant BCGs. 
Specifically, we parameterize the magnitude-redshift rela- 
tion as 

m obs = M+51og 10 (^ L /10pc)-^loglO(( I ^)(i r ^)) a + J B 

(1) 

The third term is the surface brightness correction which 
decreases the scatter by a factor of ~2. The last term pa- 
rameterizes evolution and is based on results from Bruzual 
& Chariot models (Bruzual & Chariot 1993). Because the 
magnitude changes nearly linearly with redshift, the na- 
ture of the parameterization is not critical (a linear func- 
tion of z is on average only different by 0.01 mag from 
our adopted relation). Here, and throughout this work, 
we work in the observed frame, not the rest-frame, and 
as such, any "evolution" will include both real evolution 
in the stellar populations plus K-corrections. The free pa- 
rameters in this fit are A, B, a, and 7 (see Gonzalez et 
al. 2001 for a full discussion). We fix the cosmological pa- 
rameters VLm =0.3 and J7a = 0.0 for computing Dl. This 
choice of cosmology has no impact on the results, as we are 
only interested in using this empirically calibrated relation 
to estimate redshifts within the redshift range spanned by 
the calibration data. 

This parameterization of the relationship works well. A 
comparison between zbcg vs. Zspec 

is presented in the 
lower panel of Figure ||. The rms about the 1:1 line is 
0.05. Because the errors in z spec are negligible in compar- 
ison to those in Zbcg, we adopt <J ZBCG =0.05. In addition 
to our high redshift cluster sample, there are three previ- 
ously known clusters that fall within our survey area and 
have published redshifts. We measure their BCG magni- 
tudes from the drift scans (open circles in Figure ||) and 
find that they are in remarkably good agreement with our 
parameterization of the tcibcg vs. z relation, even though 
they lie outside the redshift calibration range. For all of 
our cluster candidates with detected BCGs, we invert this 
relation to estimate a cluster redshift. Uncertainties due 
to misidentification of the BCG or on the the measurement 
of E are discussed by Gonzalez et al. (2001) who conclude 
that (Jzbcg 1S 0-08. The quoted uncertainties on the pho- 
tometric redshift indicators may be underestimates of the 
uncertainties for the catalog as a whole, because the spec- 
troscopic clusters are generally richer cluster candidates, 
which may have more homogeneous BCGs. 

3.2. Estimating Redshifts Using the Cluster Galaxy 
Luminosity Function 

The galaxy luminosity function can be characterized us- 
ing a Schechter (1976) function, with two free parameters 
m*j and a. Our photometric redshift estimation technique 
relies on \ 2 minimization of this functional form to our 

those objects are near bright stars. Subsequent catalogs (including our 
using larger masked regions around bright stars and galaxies. Because a 
there are two clusters with spectroscopic redshifts that are not include 



/-band data to measure mj . The /-band data are less 
sensitive to stellar population differences among galaxies 
than bluer optical bands, and much more observationally 
efficient than IR bands. Because our data at best reach 
only a few magnitudes below m} , we fix the slope of the 
faint end of the luminosity function at a — —1.25 (Lugger 
1986) and only use galaxies with mr <J 23 in the fit. 

To generate the cluster luminosity function, we define a 
cluster center and radius, and remove contamination due 
to projected stars and interloping galaxies. The initial 
cluster center is taken to coincide with the centroid of the 
.7 LSB fluctuation in the drift scan detection. The centroid 
of the positions of the 30 galaxies closest to the initial 
cluster center is used as the final center. The overall shift 
from initial to final center is typically less than 20 pixels 
(<; 10"). We fit the surface number density of galaxies 
within 700 pixels of the cluster center with a radial profile 
of the form C\ + Ca/(1 + r ) 2 ■ The surface number den- 
sities for our 12 spectroscopically confirmed clusters with 
/-band luminosity functions are shown in Figure ||. The 
initial cluster radius is defined to be the radius at which the 
surface density of galaxies is 10% greater than the back- 
ground surface density, where the background is defined 
at r > r c i + 100 pixels, and r c i is the cluster radius. Be- 
cause our data are primarily sensitive to the bright cluster 
galaxies that lie in the central regions of the clusters, the 
radial fits yield a rather crude estimation of the cluster 
radius. Therefore, the initial cluster/background radius 
is incrementally increased by 10% until an increase in ra- 
dius adds more background galaxies than cluster galaxies. 
We typically become dominated by background galaxies 
at r c i ~ 150 — 300ft. -1 kpc and consequently our lumi- 
nosity functions are derived from central cluster galaxies 
only. Stars are removed from the object catalog using SEx- 
tractor's star/galaxy separation index, CLASS- STAR. To 
account for the effect of differing seeing conditions on the 
classification parameter, we examine the behavior of peak 
surface brightness of objects vs. CLASS- STAR, magni- 
tude vs. CLASS- STAR, and magnitude vs. peak surface 
brightness. We find that objects with CLASS- STAR ;> 
0.95 are typically stellar-like based upon their locus in the 
magnitude vs. peak surface brightness diagram and a vi- 
sual inspection of their isophotal shapes. (CLASS- STAR 
— > 1.0 are stars and CLASS-STAR — > 0.0 are galaxies). 
Because our clusters typically cover less than 2% of the to- 
tal area of an image, we use individual images to statisti- 
cally subtract background galaxies. By using our own im- 
ages rather than published background counts, we sidestep 
such issues as differences in completeness and variations in 
counts due to projected large-scale structures. 

We vary m} and search for the minimum x 2 to find 
the best fitting Schechter function. We apply this proce- 
dure using two different cumulative luminosity distribu- 
tions. The first distribution is the cumulative number of 
galaxies as a function of magnitude, while the second is the 
cumulative luminosity of galaxies as a function of magni- 
tude. The first distribution weights the faint end more 
heavily (there are more faint galaxies and therefore nec- 
essarily more bins and more weight at the faint end) and 
is consequently less affected by statistical variations. The 

inal catalog presented by Gonzalez et al. 2001) are drawn from analyses 
e utilize a surface brightness correction in our msc'G vs - 2 relationship, 
I in the calibration. 
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second distribution, on the other hand, weights the bright 
end more heavily where the errors in our photometry are 
small. The fitting begins at a magnitude corresponding 
to the brightest galaxy remaining after background sub- 
traction (presumably the BCG). Defining the faint end 
of the fit is more difficult because each cluster has a dif- 
ferent completeness limit. Rather than determining the 
completeness of each image, we exploit the fact that the 
fits to our two luminosity functions are separately more 
sensitive to the bright and faint ends. Therefore, changes 
in the faint end limiting magnitude affect one much more 
strongly than the other. The faint end limit is initially 
set to be m/ = 23 (none of our data is complete to this 
magnitude) . This limit results in a severe underestimation 
of the redshift based on fits to the cumulative number of 
galaxies. Fits to the cumulative luminosity of galaxies also 
underestimate the redshift, but to a far lesser degree. We 
then decrease the faint end limiting magnitude until the 
two redshift estimates agree to better than 15%. The aver- 
age of these two estimates is taken to be the redshift, Zi um , 
of that cluster. In Figure |4] we plot the cumulative number 
of galaxies as a function of magnitude for our 12 spectro- 
scopically confirmed clusters with /-band data. The solid 
line is the best fitting Schechter function corresponding to 

Zium • 

We calibrate the m^ - z relation using clusters with spec- 
troscopic redshifts (upper panel of Figure |5|). We parame- 
terize the relation with a linear function, nij = Ci + Q>iZ, 
with best fitting parameters Ci = 17.40 ± 0.11 and C2 = 
3.28 ± 0.24 (the rms about the fit is 0.28). The plotted 
error bars represent a mean of the data binned in redshift 
intervals Az = 0.1, where the mean is taken to be the value 
of the fit corresponding to the center of the bin. The lower 
panel of Figure || compares z\ urn and z spec . Because the 
uncertainties in z/ um are much larger than those in z spec , 
we adopt the rms scatter about the line, 0.06, as the un- 
certainty in Zl urn . 

3.3. Color Envelope 

O'Connell (1988) defined the "red envelope" as the nar- 
row locus of red galaxies in the color-magnitude diagrams 
of intermediate redshift clusters. Subsequently the exis- 
tence of the red envelope has been confirmed for clusters 
extending to redshifts of z ~1 (A93; Stanford et al. 1995; 
Lubin 1996; SED98; Kodama et al. 1998). Smail et al. 
(1998) studied 10 x-ray luminous clusters in the redshift 
range z = 0.22 — 0.28 and found a small spread in the ob- 
served (B — I) colors (~0.04) for the brighter cluster ellip- 
ticals. These authors suggest that the location of the red 
envelope in clusters would be an efficient, albeit somewhat 
crude, indicator of redshift provided the relation could be 
calibrated with clusters spanning the redshift range of in- 
terest. 

We exploit the small scatter in the CM relation and 
use the location of the red envelope in (V — I) as a red- 
shift estimator, calibrating the relation with a subset of 
our clusters that have spectroscopic redshifts. The loca- 
tion of the red envelope is not significantly affected by 
contamination from non-cluster galaxies, so we include all 
galaxies within 350/i _1 kpc projected radius of the cluster 
center. For clusters without spectroscopy, we use our esti- 
mated redshift from the luminosity function fitting, z; um , 
to determine the projected physical radius. However, the 



resulting location of the color envelope is not sensitive to 
relatively small changes in radius, thereby providing an 
independent redshift indicator. 

Wc define the location of the color envelope to corre- 
spond to the maximum change in the number of galaxies 
as a function of color. This technique can be misled by 
the presence of a tight clump of faint, blue galaxies in 
the CM diagram. To circumvent this possibility, and be- 
cause the red envelope is dominated by bright ellipticals, 
we only consider galaxies with ni/ < 22. The red envelope 
position is determined by an automated routine and sub- 
sequently checked by a visual inspection of the cluster's 
CM diagram. Figure ^, upper panel, shows the location 
of the red envelope in (V — I) as a function of redshift 
for the subset of our clusters that have spectroscopic red- 
shifts, along with clusters from A93, S97, and SED98. We 
fit the empirical relation with a parabolic function of the 
form (V — I) =Ci+C2Z+C3Z 2 , where the best fitting pa- 
rameters are Ci = -0.56 ± 0.53, C 2 = 8.50 ± 1.82, and 
C 3 = -5.44 ± 1.50 (rms about the fit is 0.16). As for 
the mj vs. z relation, the error bars represent o~ mean of 
the data binned in redshift intervals Az = 0.1, where the 
mean is taken to be the value of the fit corresponding to 
the center of the bin. At high redshift (z ;> 0.6) there are 
only small differences in the V — I colors of bright cluster 
ellipticals and the empirical relation flattens. Therefore, 
clusters with z co \ J> 0.6 should be regarded as lower limits 
to the true cluster redshift. In the lower panel of Figure [l(] 
we plot z co i vs. z spec . The uncertainties in z spec are negli- 
gible compared to those in z co i and therefore we set tr Zcoi 
equal to the rms about the line, 0.07. However, because 
the bulk of the clusters used in this determination have 
Zspec <i 0.6 where the empirical relation is steep, this er- 
ror estimate does not apply to high redshift clusters whose 
value of z co i is merely a lower limit. The small number of 
confirmed clusters at high redshift does not allow us to es- 
timate the errors associated with the lower limiting values 
of z co i. 

Color magnitude diagrams for all candidates with 
follow-up V and / data that we consider to be true clus- 
ters are shown in Figure ^. These diagrams include all 
galaxies within 350ft. -1 kpc and are not background sub- 
tracted. The clusters are shown in order of increasing red- 
shift according to zi um , if z spec is unavailable. The solid 
line shows the location of the red envelope as determined 
by the automated procedure described above. The dashed 
line corresponds to the expected location of the red en- 
velope using our derived (V — /) vs. z relation and the 
cluster's known or independently estimated redshift. For 
the majority of clusters, the agreement between the two 
estimates is excellent. 

3.4. Classification of Cluster Candidates and Their 
Photometric Redshifts 

Without spectroscopy of a number of galaxies per cluster 
field, the confirmation of candidate clusters is indirect. We 
use the following criteria as guidelines in the classification 
process: 1) the appearance of the background subtracted 
/—band luminosity function, in particular whether there 
is a well populated and defined bright end of the luminos- 
ity function, 2) the presence of a prominent red envelope 
in the (V — I) CM diagrams, and 3) a clear peak in the 
surface density of galaxies, in particular a high degree of 
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spatial clustering of red galaxies on the image. Any cluster 
candidate that strongly satisfies any of these three criteria 
is considered a viable cluster. However, we do not require 
successful candidates to satisfy all three criteria for two 
reasons. First, we do not have the necessary data to eval- 
uate every cluster equivalently. Second, we have confirmed 
clusters that do not strongly satisfy all three criteria (for 
example, high-z clusters generally do not have prominent 
red envelopes due to the combination of cosmological shift- 
ing of light out of the y-band and our modest exposure 
times). 

Due to the heterogeneous nature of the data set and 
the intermediate stage of our cluster selection algorithms, 
we acknowledge that our classification procedure is sub- 
jective. However, we demonstrate that random fields and 
non-cluster detections (LSBs and spurious detections) do 
not yield well-defined luminosity functions nor prominent 
red envelopes. In Figure || (panels a-f) we plot luminos- 
ity functions of random fields on our images. Negative 
values are arbitrarily assigned a value of 0.04. The ran- 
dom centers are chosen to lie entirely outside the cluster 
region. The random cluster radius is determined from sur- 
face density profiles of the random fields in the same man- 
ner as for clusters. Because these fields do not contain 
sharply peaked overdensities of galaxies, their fitted radii 
are generally much larger than our typical cluster radii. 
Figure || clearly demonstrates that random fields do not 
produce sensible luminosity functions. Next we plot the 
luminosity functions of three candidates classified as either 
an LSB galaxy (panel g), a spurious detection (panel h), 
and a doubtful cluster (panel i). These three objects are 
typical of our non-cluster detections in that their luminos- 
ity functions are either very steep (leading to zi um > 1), 
noisy, or flattened (leading to zi um <^ 0.25). The third 
object is a good illustration of a borderline case (Figure |[ 
panel i). The luminosity function is too uncertain to use 
as the sole criteria for classification as a cluster. There- 
fore, because this object lacks both color information and a 
strong clustering of galaxies near the LSB detection peak, 
we exclude it from our current analysis. Deeper follow-up 
imaging of this candidate is necessary to finalize its clas- 
sification. Next, we examine the CM diagrams of random 
fields and non-cluster detections. For reference, we repro- 
duce the CM diagram of cll404— 1216, a spectroscopically 
confirmed cluster at z — 0.39 (Figure |[ panel a). We 
choose 5 random centers on the image that lie entirely 
outside the cluster region (r > 700 h^ 1 kpc) and consider 
all galaxies within 350/i _1 kpc of the random center to 
generate CM diagrams of random fields (Figure 0, panels 
b-f). The solid line shows the maximum change in the 
number of galaxies as a function of color as determined by 
our automated procedure. The random fields do not show 
any coherent structure that might be mistakenly identified 
as a red envelope. Similarly there is no indication of the 
presence of a red envelope in the CM diagrams for three 
rejected cluster candidates (Figure |, panels g-i). 

Based upon our three criteria, Table 2 lists the 65 candi- 
dates that we classified as clusters, along with z speci z\ um , 
z C oi, zbcg, m Jj the location of the red envelope in (V — I) 
and (I—K' ). Only the results using these clusters are 
presented in the photometric redshift calibrations previ- 
ously discussed. Of our 65 clusters, 58 (92%) have zi um , 



36 (57%) have zbcg, 30 (48%) have z coh and 17 (27%) 
have z spec . Figure [f(J compares the estimated redshifts 
derived using the three photometric redshift estimation 
techniques: z /um vs. z cot (left panel), z lum vs. z B cg (cen- 
ter panel), and z co i vs. zbcg (right panel). All three 
techniques correlate with >99% confidence according to 
the Spearman rank test and have rms scatters about the 
1:1 lines that are in agreement with the errors adopted for 
each photometric redshift estimator. Both of these factors 
convince us that these objects are bonafide clusters. Of 
the 112 cluster candidates for which we obtained follow- 
up imaging, we rejected 49 objects (44%). This relative 
success or failure rate is not indicative of the final cluster 
catalog, however, because we specifically tested the limits 
of our previous cluster selection algorithms by selecting 
some marginal candidates. 

4. DISCUSSION 

We now reverse roles and use our various photometric 
measures to study cluster galaxy evolution. In compar- 
ison to other recent studies with superior data on indi- 
vidual clusters (eg. S97, SED98), we cannot study any 
of our clusters in similar detail. However, we provide an 
independent set of clusters and a significant enlargement 
of the sample size with the combined dataset. For exam- 
ple, in the study of the evolution of M} , we increase the 
sample of clusters at these redshifts by 200% (see §4.2). 
Because we are searching for general trends, it is arguable 
that more clusters with less precise data may be a supe- 
rior sample than fewer clusters with more precise data. 
Clearly, the two types of studies are complementary. 

In such an endeavor, one is always uncertain whether 
the same class of object is being compared between high 
and low redshifts. We take complimentary approaches in 
this discussion to root out potential difficulties. We begin 
by examining the integrated population as characterized 
by the luminosity evolution of MJ as a function of red- 
shift. Next we focus on the galaxies that occupy the very 
narrow locus in color in the cluster CM diagram charac- 
teristic of elliptical galaxies. While we cannot guarantee 
that galaxies do not move in and out of any one class of 
galaxies (in fact we will later suggest that they do), we can 
at least describe how that class evolves. Measurements of 
well-specified populations at the very least provide clear 
targets for simulations. Finally, because our survey is bi- 
ased towards more massive clusters at high redshift, we 
search for correlations between these empirically measured 
cluster galaxy properties and cluster mass. 

4.1. Luminosity Evolution 

One way to characterize luminosity evolution in clus- 
ter galaxies is to parameterize changes in the value of Mj 
as a function of redshift. The value of mj is a by-product 
of our luminosity function photometric redshift estimation 
technique. If an independent determination of redshift is 
available, mj and z can be used to constrain evolution in 
Mj . The left panel of Figure [ll] shows the derived value of 
M} both for clusters with spectroscopic redshifts and those 
that have either zbcg or z co i. For clusters with both Zbcg 
and z co i, we use zbcg because the associated uncertainties 
are smaller and because zbcg is less likely than color to 
be directly tied to the overall evolution of the bright end of 
the galaxy luminosity function. There is a slight tendency 



Nelson et al. 



9 



for our spectroscopic clusters to have a brighter M*j than 
our non-spectroscopic sample, although we find that the 
difference is statistically insignificant. The galaxy magni- 
tudes arc corrected for Galactic extinction, but are not K- 
corrected. Although there is considerable scatter in MJ at 
any z, a correlation exists between MJ and z at the >99% 
confidence level according to a Spearman rank correlation 
test. We parameterize the full data set with a linear func- 
tion of the form (MJ +5 log ft.) = Ci + C 2 z, with best fit 
parameters of Ci = -21.74 ± 0.12 and C 2 = -0.88 ± 0.24 
(rms about the fit is 0.34). The error bars in Figure [ll] 
represent er mea „ for clusters binned in Az = 0.1, where 
the mean is taken to be the value of the fit corresponding 
to the center of the redshift bin. This parameterization is 
valid over 0.3< z <0.9, the redshift range probed by our 
sample, and for Sl = 0.2(fiA = 0). We conclude that the 
luminosity in the observed / frame of a characteristic clus- 
ter galaxy is 0.53 ± 0.28 + 5 log h mag brighter at z = 0.9 
than at z = 0.3. 

We test for possible systematic errors in a variety of 
ways. First, this result is consistent with that obtained 
using only clusters with spectroscopic redshifts (for which 
Mj is 0.59 ± 0.17 + 51og/i mag brighter over the same 
redshift range), which argues against the influence of po- 
tential problems with the photo-z measurements. Second, 
the result is consistent with that obtained using only the 
literature clusters, which again argues against the influ- 
ence of problems with either our photometry, the fitting 
of m}, or the photo-z's. Lastly, we find that our mea- 
surement of m*j is independent of the limiting magnitude 
of the observations, where we estimate the completeness 
limit of each field using the color-magnitude diagrams and 
identifying the / magnitude at which the red sequence is 
lost. There is no evidence that the magnitude limit affects 
our determination of m* (Figure |lj) . 

How much of this apparent evolution is due to K- 
corrections and how does our observed evolution in MJ 
compare to theoretical predictions? A direct comparison 
requires detailed modeling of the formation and evolution 
of the cluster galaxy population as a whole and is beyond 
the scope of this paper. We take a simpler approach and 
compare the evolution of M? to the expected luminosity 
evolution of both a spiral and an elliptical galaxy. We 
caution that due to the various complications involved in 
the interpretation of evolution of the luminosity function, 
principally that we are studying the full wide range of mor- 
phological classes present at these magnitudes, any con- 
clusions drawn must be viewed as merely suggestive. In 
particular, the error bars do not include uncertainties in 
the relative normalization of the models and assume that 
the scatter in the measurements is Gaussian about a single 
true value. 

The right panels of Figure [ll] compare our observa- 
tions of MJ vs. z to spectral synthesis model predictions 
obtained using the GISSEL96 code (Bruzual & Chariot 
1993). Because parameterizing data involves some arbi- 
trariness in the selection of the functional form, we choose 
to compare directly to the data, binned in redshift inter- 
vals of Az = 0.1. The associated error bars are o~ mean 
of the binned data. We consider two types of star forma- 
tion histories: 1) an initial 10 7 Gyr burst of star forma- 
tion (elliptical-like; upper panel), and 2) an exponentially 
declining star formation rate such that /i is the fraction 



of the galaxy's mass converted into stars and is given by 
1 — e _1 / r , with n =0.1 and r in Gyr (spiral-like; lower 
panel). All models have Salpeter initial mass functions 
for masses between 0.1 M and 100 Mq. The models are 
normalized to the data at z =0.45, which corresponds to 
the center of the redshift bin that has the greatest number 
of spectroscopically confirmed clusters. The models are 
computed in the observer's frame and thus can be directly 
compared to our observations. 

The most immediate conclusion is that no model using 
a spiral galaxy SED fits the data. No evolution and pas- 
sive evolution models predict roughly the same luminosity 
evolution for a spiral, regardless of Zf orm , and therefore 
only two representative cases are plotted. The model spi- 
ral galaxy's luminosity remains roughly constant over the 
redshift range 0.3 < z < 1 and is inconsistent with the ob- 
served brightening of . Turning to the elliptical galaxy 
evolution models, we find that the no evolution model is 
ruled out. Also, passive evolution models with Zf orm > 2 
are poor fits to the data. In these models, a MJ galaxy 
formed at a sufficiently early epoch that appreciable star 
formation ceased well before z so the galaxy is not 
evolving sufficiently over the redshift range probed by our 
sample to counter the K-correction. In both of these cases, 
actually dims with increasing redshift due to the cos- 
mological shifting of light out of our observed passband. 
Our observations of MJ agree best with passive evolution 
models with a relatively low Zf orm ~ 1.7. In these models, 
a M} galaxy underwent its burst of star formation rather 
recently, and we are simply tracking the decline in lumi- 
nosity due to the passive aging of the stellar populations. 
These models do not include mergers or accretion, so we 
only infer from this result that the luminous galaxies in 
these clusters are somewhat brighter than expected hat 
they had formed at high redshifts and since been quies- 
cent. 

In contrast to some earlier studies that focus on the red- 
dest optical bands or even extend into the near-IR wave- 
lengths (e.g. de Propris et al. 1999) to avoid the difficulty 
in modeling the rest-UV spectral energy distribution of 
galaxies, or that use different passbands for different clus- 
ter redshifts to minimize K-corrections (e.g. SED98), we 
have used primarily V and / (and K' for a smaller frac- 
tion of clusters) in the observed frame. We consider the 
use of relatively blue rest-frame optical data an important 
complement to existing data. Before comparing our con- 
clusions to those of previous studies, we discuss whether 
our approach compromises our study. The disadvantage 
of using V and / is that it places greater demands on the 
models because the rest-frame UV must be modeled well 
(observed V corresponds to rest-frame 2900A at z = 0.9). 
The advantage is that it probes the spectral range that is 
most sensitive to recent star formation. The models may 
be more robust at redder wavelengths, but the effect of re- 
cent star formation is more subtle. If our results agree with 
previous studies, then we can conclude that the theoreti- 
cal models are working in the [/-band and that the previ- 
ous conclusions regarding the evolution of cluster galaxies 
hold. If our results disagree with previous studies, then 
we can conclude either 1) that the models may have prob- 
lems in the [/-band, 2) that our clusters are different than 
those studied previously, 3) that evolution beyond simple 
passive evolution is fairly subtle and only detectable at 



10 



Cluster Galaxy Evolution and the LCDCS 



bluer wavelengths, or 4) that there are systematic errors 
in either our data or previous data (or both). 

Our findings are in qualitative agreement with previous 
studies of luminosity evolution. Using WFPC2 images of 
C10939 (z — 0.41) and ground-based i?-band images of 
Coma, Barrientos et al. (1996) measured the luminosity- 
effective radius relation for bright ellipticals and deter- 
mined that passive evolution models with Zf orm ~1.2 pro- 
vided the best fit to their data. S97 measured composite 
luminosity functions for early and late-type galaxies in 10 
clusters at 0.37 < z < 0.56 from WFPC2 images. The el- 
lipticals show a modest, but significant, brightening of My 
consistent with passive evolution and Zf orm ~ 3. In con- 
trast, the late type galaxies show no evidence for a bright- 
ening in My over the same redshift range. Although our 
luminosity functions consist of the entire range of mor- 
phological classes, we are limited to the bright galaxies, 
which are presumably ellipticals. While S97's luminosity 
evolution best matches passive evolution with a slightly 
higher Zf orrn than our observations, their sample only 
probes out to z =0.56. The upper right panel of Fig- 
ure [ll] demonstrates that significant differences in passive 
evolution models do not appear until z > 0.6 and thus we 
conclude that we are not in disagreement with S97. Addi- 
tionally, we are in agreement with de Propris et al. (1999) 
who find that M|^ departs from no-evolution model predic- 
tions at z > 0.4 and that this departure is consistent with 
passive evolution with intermediate Zf orm . Again, because 
of many unquantifiable uncertainties we do not highlight 
minor discrepancies between the studies, rather we search 
for broad qualitative agreements or disagreements. 

4.2. Red Envelope Evolution 

The previous approach analyzes the evolution in the 
combination of galaxy populations, therefore the results 
are somewhat difficult to interpret. Alternatively, we can 
select a single galaxy population and study its evolution. 
We investigate a group of galaxies that inhabit a well- 
defined locus in the cluster CM diagram, the red envelope. 
We plot the location of the red envelope in V — I versus 
redshift (Figure [l3] left panel) and in I—K' (Figure 14|, 
left panel). Again, we use clusters with independent pio- 
tometric redshifts (in this case zbcg or zi um ) and clusters 
with spectroscopic redshifts to parameterize the relations. 
For V—I we fit the observed relation with a parabolic func- 
tion of the form (V — I) =Ci+C2^+C3Z 2 , where the best 
fitting parameters are Ci = -0.24±0.28, C 2 = 7.42±1.03, 
and C 3 = -4.61 ± 0.91 (the rms about the fit is 0.16). 
For I—K' the observations are best fit with a quadratic 
function of the form (I-K 1 ) =Ci+C 2 z+C 3 z 2 +C 4: z 3 , 
where the best fitting parameters are Ci = 0.66 ± 0.65, 
C 2 = 9.50±3.72, C 3 = -14.72±7.01 and C 4 = 8.72±4.29 
(the rms about the fit is 0.18). Both parameterizations are 
valid over 0.3 < z < 0.9, the redshift range of our clusters. 
The error bars are a mean for clusters binned in Az = 0.1, 
where the mean is taken to be the value of the fit corre- 
sponding to the center of the redshift bin. We summarize 
our empirically derived luminosity and color evolution pa- 
rameterizations in Table 3. 

The brightest, reddest galaxies, which are presumably 
the cluster's oldest ellipticals populate the red envelope. 
We compare our measured color evolution in (V — I) (Fig- 
ure |l3|, right panel) and (I—K' ) (Figure |l4|, right panel) 



with model predictions obtained using the GISSEL96 code 
(Bruzual & Chariot 1993) for an elliptical galaxy undergo- 
ing passive evolution or no evolution. The parameters of 
the models are the same as those given above. When com- 
paring colors, in contrast to the magnitude comparisons 
of the previous section, no normalization is necessary. Be- 
cause there is a certain amount of arbitrariness in param- 
eterizing the data with a function, here we simply bin the 
data in intervals of Az = 0.1. The error bars are a mean of 
the binned data. Turning first to the evolution in (V — I), 
at low redshift the data lie intermediate between no evo- 
lution and passive evolution models with a high formation 
redshift, Zf orm >5. However, at high redshift where the 
differences in model predictions increases, we find that the 
observations are inconsistent with the no evolution predic- 
tion. Our data are most consistent with models of elliptical 
galaxies with passive evolution and Zf orm J>5. Examining 
the color of the red envelope in I—K' , we find that no 
evolution models are completely ruled out. Such models 
predict colors that are too red compared to the location 
of the red envelope at all z. These data are most consis- 
tent with a passive evolution model with a high epoch of 
formation, z form ;> 5. 

Similar conclusions have been reached by previous 
authors, using both optical and infrared-optical colors. 
RS measured galaxy colors in the rest-frame Stromgren 
u v b y filters and found a degree of blueing of their sam- 
ple from z ~ 0.4 to 0.9 that implies a redshift of formation 
of - 10. From the PDCS using the WFPC2 bands V 4 
(F555W) and J 4 (F785LP), Lubin et al. (1996) find the 
red envelope of their clusters is bluer in V4 — I4 than no 
evolution predictions, and implying passive evolution with 
Zform > 2. Using optical-infrared colors, both A93 and 
SED98 observe the blueing trend compared to no evolu- 
tion model predictions. Their observations agree best with 
models in which galaxies form at z ~ 2 to 4 and passively 
evolve thereafter. We conclude that one can identify clus- 
ter galaxies that appear to have formed the bulk of their 
stars at z ~ 5 by selecting galaxies on the red envelope. 

4.3. Cluster Selection Biases 

One concern that cannot yet be entirely addressed is 
that the comparison across a large redshift range incor- 
porates clusters with different properties. We expect to 
detect only the richest, most massive systems at high red- 
shifts and a wider range of systems at low redshifts. Be- 
cause we do not have x-ray observations of our clusters, 
we cannot yet address this potential bias directly. How- 
ever, an indirect cluster mass estimation technique indi- 
cates that indeed we only detect the most massive sys- 
tems at high redshift (see Gonzalez et al. 2001). Using 
17 known clusters with x-ray data for which we have drift 
scans, we find a correlation between the x-ray luminosity, 
h x , and the peak surface brightness of our cluster detec- 
tion, S. Using this correlation, we find that we detect 
clusters with a wide range in mass (groups to rich clus- 
ters) at low to intermediate redshift. However, at high 
redshift (z ;> 0.7) we detect only the most massive sys- 
tems (e.g. we can detect MS 1054-03 with L x ~ 9 x 10 44 
ergs s" 1 ). 

Whether the galaxy properties of these systems depends 
on the global characteristics of the cluster is not well de- 
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termined. Recent work suggests that the color evolution of 
bright ellipticals and the evolution of M* are not strongly 
affected by cluster mass. Stanford et al. (1995) studied 
an inhomogeneous sample of 19 clusters with < z < 0.9 
drawn from a variety of optical, x-ray, and radio-selected 
surveys that span a wide range in cluster mass (3.0 x 10 43 
ergs s~ 4 <L X < 2.0 x 10 45 ergs s _1 ). They found the 
color evolution of the bright elliptical galaxies in these 
clusters does not depend strongly on either L x or opti- 
cal richness. Similarly, Smail et al. (1998) found <2% 
scatter in the rest-frame UV-optical colors of the bright 
elliptical sequence in 10 x-ray luminous ROSAT clusters 
at 0.22 < z < 0.28. 

The evolution of M* also does not seem dependent upon 
cluster mass. Using a heterogeneous sample of 38 clusters 
at 0.1 < z < 1, de Propris et al. (1999) fit Schechter func- 
tions with a fixed a = —0.9 to their K-band data. They 
found no statistically significant difference in the measured 
value of m^- as a function of z when they divided their 
sample based upon x-ray luminosity or optical richness. A 
similar result is found in optical bands. Garilli, Maccagni, 
& Andreon (1999) studied a sample of 65 Abell and x-ray 
selected clusters at 0.05 < z < 0.25 in the Gunn g, r, and 
i bands. They constructed composite luminosity functions 
and found similar values of M* in all three bands in both 
the Abell and x-ray selected samples and also between 
subsamples divided by optical richness. 

We use our empirically derived correlation between the 
peak surface brightness of the cluster detection, S, and 
x-ray luminosity, L^, and search for trends between clus- 
ter galaxy properties and cluster mass. In Figure we 
compare E(l + z) 4 to (Mj +5 log ft.) (left panel) and the lo- 
cation of the red envelope in (V — I) (right panel). Neither 
of these quantities are significantly correlated according to 
the Spearman rank test. Although current work suggests 
that certain cluster galaxy properties are not strongly de- 
pendent upon cluster mass, this issue merits further study 
and requires large cluster samples such as that developed 
here. Better mass estimates for these clusters are critical 
in addressing this issue. 

4.4. Cluster Galaxy Evolution 

The various empirical results presented above, in con- 
junction with published results, present a rather complex 
picture of cluster galaxy evolution. First, the galaxies on 
the cluster E/S0 sequence appear to be evolving as ellip- 
tical galaxies with very large formation redshifts (j> 5). 
Second, the bulk of the luminous cluster galaxies (as de- 
scribed by the luminosity function) appear to be evolving 
similarly to elliptical galaxies with low formation redshifts 
(<; 2). Are these two independent galaxy populations, or 
is there one process that explains all of them? 

It is evident that there are some cluster galaxies, those 
on the red E/S0 sequence, that consist of very old stellar 
populations. The comparison of spectral synthesis models 
with the V — I and I—K' data suggests that these galaxies 
have a dominant population that formed at z > 5. De- 
tailed morphological and spectroscopic studies of clusters 
at high redshift confirm that these galaxies are a subset of 
the clusters' elliptical population (e.g. van Dokkum et al. 
1999, SED98). However, because of the non-uniqueness of 
the spectral synthesis model, it is important to differenti- 
ate between a stellar population that is predominantly old 



and one that is entirely old. Although models with recent 
formation redshifts (i.e. galaxies with stellar populations 
that arc entirely formed at z < 2) do not fit the color en- 
velope data, it is not evident from the previous discussion 
whether or not models with some recent star formation 
superposed on an older population fit the data. 

The observed behavior of the luminosity function sug- 
gests significant evolution in the cluster galaxy popula- 
tion. Because of our relatively bright magnitude limit, 
this result pertains almost entirely to luminous galaxies 
(L > L„) that primarily become (or have always been) el- 
liptical galaxies in the cores of these rich clusters. There- 
fore, on the one hand we have a seemingly quiescent (al- 
though not necessarily quiescent) elliptical population and 
on the other an evolving population. Is it possible that the 
entire population is evolving? One solution is suggested 
by the observation of red mergers in MS1054— 03 by van 
Dokkum et al. (1999). These galaxies will remain on the 
red sequence after the merger but M* will be affected. 

Can mergers that do induce some star formation also ac- 
count for the observations? In Figure |l6| we compare the 
observed evolution of the red envelope in V — I (left panel) 
and I—K' (right panel), and spectral synthesis models 
(Bruzual & Chariot 1993) of starburst episodes in ellipti- 
cal galaxies. These models consist of a passively evolving 
elliptical galaxy that formed at Zf orm = 10 and subse- 
quently underwent a burst of star formation at low red- 
shift, Zb urst — 1.5 (solid lines) and z\, urst — 1 (dashed 
lines). The duration of the burst is 10 7 years and has 
strengths (amount of galaxy's mass that is converted into 
stars) of 10% (thick lines) and 25% (thin lines). At the 
peak of the starburst, the bluening of the galaxy's inte- 
grated light is dramatic (~2.5 magnitudes in V — I), but 
also short-lived. The bulk of the post-starburst period 
involves a moderate blueing of the galaxy (~ 0.2-1 mag- 
nitude in V — I). The strength of the burst does not sig- 
nificantly affect the amount of blueing, but does affect its 
duration. It is evident that an elliptical galaxy returns 
fairly quickly to its standard location on the red E/SO 
sequence and that the galaxy appears similar to one that 
simply had a single burst of star formation at z = 10. 
We stress that as galaxies evolution they may or may not 
be included in subpopulations like the red sequence, which 
leads to complicated evolutionary effects (see van Dokkum 
and Franx 2001 for a more complete discussion of such 
matters). It is generally thought that galaxies primarily 
move onto the red sequence from a bluer part of the color- 
magnitude space (cf. Kodama & Bower 2001 for a recent 
treatment), we simply note here that the flow may go both 
ways. 

Detailed spectroscopic observations and analysis of one 
high-redshift cluster (MS1054-03 at z = 0.83; van Dokkum 
et al. 1999) provides observational confirmation of this hy- 
pothesis. They demonstrate that 17% of the L ;> L* cluster 
population is undergoing a major merger. Most of these 
galaxies will probably evolve into luminous (^2L*) ellip- 
tical galaxies. Assuming that the galaxy population of 
MS1054-03 is representative of clusters at that redshift, 
van Dokkum et al. (1999) estimate that ~50% of present- 
day cluster ellipticals were assembled in mergers since 
z ~ 1. Most of the mergers involve red spheroidal galaxies 
with no detected [Oil] A3727 emission, but some do ex- 
hibit enhanced Balmer absorption indicative of a modest 
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recent starburst. They measure the scatter in the rest- 
frame U — V colors of morphologically selected ellipticals 
and SOs to be 0.027 ± 0.013, while that of the combined 
sample of ellipticals and mergers is significantly higher, 
0.054 ± 0.011. Evolving the CM relation, the scatter of 
the combined sample drops to ~ 0.035 at z = 0.5 and 
~ 0.015 at z = which is consistent with the small scatter 
observed in the CM relation at low redshift (Bower, Lucey, 
& Ellis 1992; Ellis et al. 1997; Stanford et al. 1998; van 
Dokkum et al. 1998). Using the colors of the mergers they 
estimated their luminosity-weighted mean redshift of for- 
mation is Zf orm ;> 1.7, which is in remarkable agreement 
with our Zf orm derived from the evolution of M^ . 

The triggering of star formation in otherwise red E/SO 
galaxies has a few consequences: 1) it will weaken or 
destroy the distinctiveness of the E/SO sequence in the 
CM diagram and 2) it will systematically brighten the 
population of ellipticals throughout the duration of the 
burst. Both of these trends are observed in our data, 
and in previous observations. Morphologically, the blue 
galaxies appear to be predominantly disky (Couch et al. 
1994, 1998; Dressier et al. 1994, 1999; Oemler, Dressier, 
& Butcher 1997; Poggianti et al. 1999), but whether 
these are spirals evolving to SO's and/or SO's evolv- 
ing to E's is not addressed by our data. A further 
complication, whether these evolutionary processes are 
driven/influenced/ terminated by the cluster environment, 
is again not addressable with our data because we have no 
proper field comparison sample. 

Alternatively, one could hypothesize that a significant 
fraction of the galaxies on current-day E/SO cluster se- 
quences came not from previous early-type galaxies, but 
from some major morphological evolution of later types 
(i.e. spiral-spiral mergers). However, the tightness of the 
red sequence (cf. Smail et al. 1998) and the detailed anal- 
ysis of galaxy spectra from a sample of nine clusters at 
0.37 < z < 0.56 (Jones, Smail, & Couch 2000) demon- 
strate that these galaxies could not have experienced a 
major episode of star formation in their relatively recent 
past. The latter study concludes that the progenitors of 
early-type cluster galaxies are themselves early-type clus- 
ter galaxies. 

In closing, our results do not drastically alter trends seen 
previously from smaller cluster samples. Instead, they pro- 
vide better sampled data with which to measure evolution- 
ary trends in cluster galaxies for redshifts from 0.3 to 0.9. 
The large sample provides confirmation that these trends 
are not dominated by a one or two "odd" clusters at high 
redshift. The empirical results are straightforward (albeit 
open to a range of interpretations) and provide tests of 
theoretical models. 

5. SUMMARY 

We analyze photometry and spectroscopy of a sample 
of 63 clusters at 0.3 < z < 0.9 drawn from the Las Cam- 
panas Distant Cluster Survey. Seventeen of our clusters 
have spectroscopically confirmed redshifts and the remain- 
ing 46 have photometrically estimated redshifts. Using 
deep V , I, and K' follow-up imaging we measure the lu- 
minosity and color evolution of cluster galaxies in order to 
empirically constrain models of cluster galaxy formation 
and evolution. Our primary results are: 
1) We develop three photometric redshift estimation tech- 



niques based upon the magnitude of the brightest cluster 
galaxy measured directly from our drift scan images, the 
7- band luminosity function, and the location of the red 
envelope in both V — I and I— K' . All three techniques 
work well with internally estimated errors of cr ZBCG = 0.05, 
o'lum — 0.06, and <j co i = 0.07. These uncertainties are de- 
rived from the spectroscopic sample and may not represent 
the uncertainties for the entire catalog. 

2) We find no systematic difference in the photometric 
properties (i.e. luminosity function, location of the red 
envelope, fraction of blue galaxies) of our spectroscopi- 
cally confirmed clusters and the rest of the sample. Nor 
do we find any systematic differences between the proper- 
ties of our clusters and those published in the literature. 
Both of these factors lead us to believe that our cluster 
finding algorithm and our cluster classification criteria are 
successful. 

3) Using 44 of our clusters and 19 clusters from the lit- 
erature (A93; S97; SED98) for which /-band data are 
available, we parameterize the redshift dependence of Mj 
in the observed frame: M* = (-21.74 ± 0.12) - (0.88 ± 
0.24)z - 5 log h for 0.3 < z < 0.9, 

4) By combining the V and / photometry of 30 of our 
clusters and 14 clusters from the literature (A93; S97; 
SED98) we parameterize the redshift dependence of the 
V— I color of the E/SO red sequence in the observed frame: 
[V -I) = (-0.24±0.28) + (7.42±1.03)z-(4.61±0.91)z 2 
for 0.3 < z < 0.9, 

5) By combining 13 of our clusters with 15 clusters from 
the literature (A93; SED98) for which / and K' data 
are available, we parameterize the redshift dependence of 
the I—K' color of the E/SO red sequence in the observed 
frame: {I-K 1 ) = (0.66 ± 0.65) + (9.50 ± 3.72)z - (14.72 ± 
7.01)z 2 + (8.72±4.29)z 3 

6) Using the peak surface brightness of the cluster detec- 
tion fluctuation, S, as a proxy for cluster L x , we find no 
correlation between E and (Mj +51og/i) or the location 
of the red envelope in (V — I). 

7) We speculate that all of these observations can be ex- 
plained with a model in which luminous early type galaxies 
(or the progenitors of current day early type galaxies) form 
the bulk of their stellar populations at high redshifts (j> 5) 
and in which many of these galaxies, if not all, experience 
a short term episode of star formation at lower redshifts 
(1.5 < z < 2). 

This sample represents a significant increase in the num- 
ber of known high redshift clusters. Although our data are 
too crude to verify if our proposed evolutionary scenario is 
correct, the statistical nature of this study provides con- 
vincing evidence that the observed trends are universal in 
massive clusters, not the result of a few outliers. The next 
step is to utilize the large sample to select several repre- 
sentative clusters of different masses and environments at 
various epochs to investigate in detail the morphologies 
and star formation rates of their member galaxies. 
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Fig. 1. — The redshift distribution of galaxies in all cluster fields observed with LRIS on Keck. Each panel represents the results for one 
candidate cluster field. Each circle represents one galaxy redshift. Filled circles represent galaxies within 1000 km sec -1 of another galaxy. 
Circles are stacked vertically within each panel to avoid overlap in the Figure. The two empty panels do not represent candidate cluster fields. 
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spec 



Fig. 2. — upper panel: The derived BCG magnitude observed in the W filter, mw , vs. z spec for clusters from our survey (filled circles) 
and three known low redshift clusters that fall in our cluster survey area and have redshifts published in the literature (open circles). Error 
bars represent the photometric errors and errors due to the surface brightness correction term, summed in quadrature. The solid line is our 
parameterization of this relation (see text for details), lower panel: Comparison of z apec vs. zbcg- The line is the expected 1:1 correlation 
of Zspec and zbcg- The data have an rms dispersion of 0.05. 
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Fig. 3. — The surface number density profiles for our spcctroscopically confirmed clusters with /-band luminosity functions. All galaxies 
within 700 pixels of the cluster center are considered. We fit the profiles with a radial profile of the form C\ + Cij{\ + r) 2 (solid lines). 
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Fig. 4. — The cumulative number of galaxies as a function of magnitude for our spectroscopically confirmed clusters with /-band data. 
Overlaid is the best fitting Schechter function (solid line). 
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Fig. 5. — upper panel: The derived value of vs. z spec for clusters from our survey, A93, S97, and SED98. The line is a linear fit to the 
data. The error bars are cr mea n for clusters binned in Az = 0.1, where the mean is taken to be the value of the fit corresponding to the center 
of the redshift bin. lower panel: Comparison of zi um vs. z apec . The line is the expected 1:1 correlation between the two redshifts indicators 
and has an rms dispersion of 0.06. 
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Fig. 6. — upper panel: The location of the red envelope in (V — I) vs. z sp ec for clusters from our survey, A93, S97, and SED98. Error bars 
are <Tmea« for clusters binned in Az = 0.1, where the mean is taken to be the value of the fit corresponding to the center of the rcdshift bin. 
lower panel: Comparison of z spec vs. z co ;. The line is the expected correlation of z spec and z co ; and has an rms dispersion of 0.07. 
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Fig. 7. — Optical color magnitude diagrams for our clusters, ordered in increasing redshift (z; um is used if z spec does not exist). Solid lines 
are the location of the red envelope as determined by our automated procedure (see text for details). The dashed line corresponds to the 
expected location of the red envelope based on the cluster's redshift and our (V — I) vs. z relation. 
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Fig. 7. — continued. 
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Fig. 8. — Cumulative luminosity functions for random fields on our images (panels a-f) and non-cluster detections (panels g-i). See text 
for details. 
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Fig. 9. — Color-magnitude diagrams for random fields on our images (panels b-f) and non-cluster detections (panels g-i). Wc reproduce the 
CM diagram of cluster 1404—1216, a spectroscopically confirmed cluster at z = 0.39, for reference (panel a). The location of the maximum 
change in the number of galaxies as a function of rcdshift as determined by our automated routine is indicated by solid lines. 
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Fig. 10. — Comparison of the estimated redshifts derived using the three photometric redshift estimation techniques, z; um , z co i, and zbcg- 



Nelson et al. 



25 



-21 



-22 - 



C30 

o 

LO 

+ 



■23 - 



-24 



~i — I — i — i — r 



A 



▲ 
▲ 



* A as a 




• z 



spec 



A Z 



phot 



oSED98 

*S97 

□A93 



-21 



-22 



o 

I — I 
LO 



. - -23 
2 



-24 
-21 



o 

lO 
+ 



-22 - 



- -23 



-24 



-| — i — I — i — i — i — I — i — i — i — I — i — i — r 



elliptical-like 







r 



p.e 

p.e. z 

p.e 



form 


1.5 


form - 


1.7 


form - 


2 



no evolution 



spiral-like 



— I— 



P e - Z f rm= 15 

no evolution 



0.4 



0.6 
z 



0.8 



0.4 



0.6 
z 



0.8 



Fig. 11. — left panel: The redshift dependence of M| vs. z (in the observed / frame) for our clusters with spectroscopic and photometric 
rcdshifts. Also included are clusters from A93, S97, and SED98. Error bars are a-mean for clusters binned in Az = 0.1, where the mean is 
taken to be the value of the fit corresponding to the center of the redshift bin. right panel: Comparison of our observed values of M* with the 
expected luminosity evolution of both an elliptical (upper) and spiral (lower) galaxy using spectral synthesis models from Bruzual & Chariot 
(1993). The data are binned in Az =0.1 and have errors bars which are cr mean of the binned data. Both no evolution and passive evolution 
models are considered. The models are normalized at z =0.45 which corresponds to the center of the redshift bin with the greatest number 
of spectroscopically confirmed clusters. 
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Fig. 12. — Comparison between the observational limiting magnitude and the measured m*. 
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Fig. 13. — left panel: The location of the red envelope in V — I vs. 2 for our clusters with spectroscopic and photometric redshifts. Also 
included are clusters from A93, S97, and SED98. Error bars arc a m ean for clusters binned in Az = 0.1, where the mean is taken to be the 
value of the fit corresponding to the center of the redshift bin. right panel: Comparison of the observed location of the red envelope in V — I 
with the expected color evolution of an elliptical galaxy using spectral synthesis models from Bruzual & Chariot (1993). The data are binned 
in Az = 0.1 and have errors bars which are cr mea n of the binned data. Both no evolution and passive evolution models are considered. 
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Fig. 14. — left panel: The location of the red envelope in I—K' vs. z for our clusters with spectroscopic and photometric redshifts. Also 
included are clusters from A93 and SED98. Error bars are <x mean for clusters binned in /\z — 0.1, where the mean is taken to be the value 
of the fit corresponding to the center of the redshift bin. right panel: Comparison of the observed location of the red envelope in I—K' with 
the expected color evolution of an elliptical galaxy using spectral synthesis models from Bruzual & Chariot (1993). The data are binned in 
Az = 0.1 and have errors bars which are cr me an of the binned data. Both no evolution and passive evolution models are considered. 
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Fig. 15. — The peak surface brightness of the cluster detection, S(l + z) 4 , compared to (M| +51ogfc) (left panel) and the location of the 
red envelope in (V — /) (right panel). The values of L x are inferred from our empirically determined correlation between S(l + z) 4 and L x ; 
they are not directly observed. Neither of the galaxy properties are significantly correlated with S(l + z) 4 according to the Spearman rank 
test. 
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Fig. 16. — The observed location of the red envelope in V — I (left panel) and I—K' (right panel) compared to spectral synthesis models 
of Bruzual & Chariot (1993) for a passively evolving elliptical galaxy that includes a secondary burst of star formation. The data are binned 
in redshift intervals of Az = 0.1 and have errors bars which are cr mean of the binned data. In all models, the elliptical forms at Zf orm = 10 
and passively evolves until experiencing a secondary burst of star formation at z^ urst = 1.5 (solid lines) or Z(, urs t = 1 (dashed lines). We 
investigate burst strengths (the percentage of the galaxy's mass that is converted into stars) of 10% (thick lines) and 25% (thin lines). 
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Table 1 
Observations 





RA 


DEC 


V 


Seeing 


I 


Seeing 


K' 


Seeing 


Cluster 


(JD2000) 


(JD2000) 


(min) 


\ j 


(min) 


c) 

V / 


(min) 


c) 


0915+4738 


09:15:51.94 


+47:38:19.9 


50 2 


1.8 


33 2 


1.9 


56 2 


2.0 


0936+4620 


09:36:06.28 


+46:20:43.2 


80 2 


2.5 


80 2 


1.9 


64 2 


2.0 


0944+4732 


09:44:21.83 


+47:32:42.6 










40 4 


2.0 


1002-1245 


10:02:01.47 


-12:45:35.1 


120 1 


1.4 


60 1 


1.4 






1002-1247 


10:02:27.14 


-12:47:13.1 


120 1 


1.4 


60 1 


1.4 






1005-1147 


10:05:43.60 


-11:47:43.1 






60 1 


1.2 






1005-1209 


10:05:49.72 


-12:09:36.5 






45 1 


1.4 






1006-1222 


10:06:29.25 


-12:22:13.7 






45 1 


1.4 






1006-1258 


10:06:18.79 


-12:58:12.5 






68 1 


1.1 






1007-1208 


10:07:42.60 


-12:08:36.0 






40 1 


1.3 






1012-1243 


10:12:14.58 


-12:43:10.9 






45 3 


0.8 


32 3 


0.8 


1012-1245 


10:12:44.35 


-12:45:37.9 


30 3 


0.9 


30 3 


0.7 






1014-1143 


10:14:56.31 


-11:43:08.1 






45 1 


1.2 






1015-1132 


10:15:19.47 


-11:32:55.5 






60 1 


1.3 






1017-1128 


10:17:45.31 


-11:28:07.8 






60 3 


1.1 






1018-1211 


10:18:46.45 


-12:11:52.8 


40 3 


1.2 


40 3 


0.6 






1023-1303 


10:23:10.11 


-13:03:52.0 






30 3 


0.6 






1024-1239 


10:24:44.94 


-12:39:55.6 


40 3 


1.1 


40 3 


0.9 






1025-1236 


10:25:08.83 


-12:36:20.0 


90 1 


1.5 


115 1 


2.7 






1027-1159 


10:27:26.31 


-11:59:33.6 


90 1 


1.5 


105 1 


1.8 






1031-1244 


10:31:50.26 


-12:44:27.2 


40 3 


1.0 


60 3 


0.8 






1032-1229 


10:32:04.91 


-12:29:43.8 


120 1 


1.8 


115 1 
15 3 


1.8 
1.0 






1041+4626 


10:41:03.79 


+46:26:36.3 


40 2 


2.5 


60 2 


2.2 


118 2 


2.0 


1059+4737 


10:59:38.03 


+47:37:38.6 


50 2 


3.3 


50 2 


2.0 


55 2 


2.0 


1100+4620 


11:00:57.36 


+46:20:38.3 










88 2 


2.0 


1136-1136 


11:36:31.95 


-11:36:07.7 


30 3 


1.3 


30 3 
20 1 


1.1 
1.1 






1136-1145 


11:36:47.07 


-11:45:33.3 


13 3 
40 3 


1.1 
1.0 


20 3 
40 3 


0.9 
0.8 


44 3 


0.8 


1136-1252 


11:36:33.54 


-12:52:03.3 


20 3 


0.8 


20 3 
20 3 


0.6 

1.4 


20 4 


2.0 


1138-1142 


11:38:06.59 


-11:42:10.3 


20 3 
60 3 


1.1 
1.3 


30 3 
60 3 


0.7 

1.2 


32 3 


0.8 


1138-1225 


11:38:14.24 


-12:25:53.9 


90 1 


1.2 


120 1 


2.0 






1100 1000 

1138 — 12z8 


1 1 . O . A A 

11:38:44.70 


— 1/:28:25.4 


30 


0.7 


o a3 
zU 


O.o 






1139-1154 


11:39:30.47 


-11:54:25.0 


30 3 


0.9 


30 3 


0.8 


10 4 


2.0 


1139-1217 


11:39:56.84 


-12:17:19.8 


20 3 


0.7 


20 3 


0.8 


32 3 


0.8 


1145-1155 


11:45:22.36 


-11:55:52.1 


20 3 
40 3 


0.7 
1.2 


20 3 
30 3 


1.0 
0.8 


50 4 


2.0 


1147-1252 


11:47:16.98 


-12:52:04.7 


20 3 
40 3 


1.0 
1.2 


20 3 
30 3 


0.9 
0.7 


32 3 


0.8 


1149-1159 


11:49:22.14 


-11:59:19.3 


3 3 


0.9 


10 3 


0.9 






1149-1246 


11:49:06.18 


-12:46:06.3 


20 3 
40 3 


0.9 
1.0 


20 3 
40 3 


0.8 
0.7 


32 3 


0.8 


1208-1151 


12:08:26.67 


-11:51:25.6 


90 1 


1.4 


68 1 


1.2 






1210-1219 


12:10:12.73 


-12:19:06.9 


60 1 


1.4 


60 1 


1.4 
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Table 1 
Observations: Continued 





R A 


T)F,C 


V 




T 

_L 


Sf Ti O' 

Ov.A.111^, 


K' 


TI (T 


V-.'l LAO L'L^A 


('TD2000 s l 


nD2000 s l 


(min) 


\ I 




\ ) 




(") 


1211-1220 


12:11:04.16 


-12:20:47.7 






45 3 


0.8 






1215-1252 


12:15:41.08 


-12:52:59.7 






68 1 


1.1 






1216-1201 


12:16:45.10 


-12:01:17.3 






45 1 


1.3 






1219-1154 


12:19:34.88 


-11:54:22.9 






60 1 


1.3 






1219-1201 


12:19:44.49 


-12:01:35.5 






60 1 


1.3 






1221-1206 


12:21:46.20 


-12:06:12.7 






45 1 


1.2 






1230+4621 


12:30:16.26 


+46:21:17.1 


67 2 


3.8 


33 2 


2.2 


72 2 


2.0 


1326-1218 


13:26:12.66 


-12:18:22.5 


40 3 


0.9 


50 3 


0.9 


36 3 


0.8 


1327-1217 


13:27:57.33 


-12:17:16.7 


20 3 
40 3 


1.1 
1.3 


20 3 
50 3 


1.0 
0.8 


42 3 


0.8 


1329-1256 


13:29:11.37 


-12:56:22.0 


90 1 

45 3 


1.6 
1.1 


95 1 
30 3 


1.4 
1.0 


20 4 


2.0 


1333-1237 


13:33:01.92 


-12:37:17.0 


20 3 


0.9 


30 3 


1.2 






1404-1216 


14:04:47.20 


-12:16:21.4 


90 1 


1.5 


60 1 


1.3 






1405-1147 


14:05:11.38 


-11:47:08.6 






60 1 


1.4 






1406-1232 


14:06:36.54 


-12:32:39.7 






45 1 


1.4 






1408-1209 


14:08:17.86 


-12:09:27.0 






68 1 


1.0 






1408-1216 


14:08:45.95 


-12:16:08.8 






68 1 


1.0 






1408-1218 


14:08:50.62 


-12:18:14.8 






68 1 


1.0 






1412-1150 


14:12:32.56 


-11:50:16.2 






8 1 
30 3 


1.2 
0.8 






1412-1222 


14:12:25.92 


-12:22:53.3 






8 1 

45 3 


1.3 
0.8 






1412-1222.1 


14:12:31.26 


-12:22:15.8 






45 3 


0.8 






1413-1244 


14:13:08.59 


-12:44:11.6 






40 1 


1.2 






1416-1143 


14:16:45.13 


-11:43:40.0 






40 1 


1.3 






1422+4622 


14:22:24.18 


+46:22:39.7 


80 2 


1.7 


83 2 


2.4 


35 4 


2.0 


1519+4622 


15:19:54.27 


+46:22:20.4 






15 2 


2.3 


45 4 


2.0 



^Las Campanas lm 
( 2 )Palomar 1.5m 
( 3 ^Las Campanas 2.5m 
W Lick 3m 
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Table 2 
Redshift Data 



Cluster 


z spec 




ov, 


z coZ 




z bcg 




mj 




V-I 




I- 


K 


&I-K 


0915+4738 


0.40 


0.43 


0.06 


0.50 


0.07 






18.77 


0.08 


2.12 


0.05 


2 


88 


0.03 


0936+4620 


0.54 


0.49 


0.06 










18.97 


0.08 






2 


96 


0.03 


0944+4732 


0.58 




























1002-1245 


0.52 


0.64 


0.06 


0.50 


0.07 


0.43 


0.05 


19.65 


0.11 


2.32 


0.03 








1002-1247 




0.72 


0.06 


0.55 


0.07 


0.57 


0.05 


19.82 


0.14 


2.48 


0.03 








1005-1147 




0.89 


0.06 






0.53 


0.05 


20.24 


0.21 












1005-1209 




0.45 


0.06 










18.82 


0.08 












1006-1222 




1.01 


0.06 






0.72 


0.05 


20.54 


0.21 












1006-1258 




0.73 


0.06 






0.47 


0.05 


19.88 


0.14 












1007-1208 




0.49 


0.06 






0.41 


0.05 


18.99 


0.08 












1012-1243 




0.49 


0.06 






0.44 


0.05 


18.99 


0.08 












1012-1245 




0.56 


0.06 


0.59 


0.07 






19.28 


0.07 


2.56 


0.07 








1014-1143 




0.30 


0.06 










17.95 


0.12 












1015-1132 




0.60 


0.06 






0.45 


0.05 


19.43 


0.07 












1017-1128 




0.50 


0.06 






0.66 


0.05 


19.04 


0.08 












1018-1211 




0.48 


0.06 


0.41 


0.07 


0.47 


0.05 


18.94 


0.08 


2.00 


0.07 








1023-1303 




0.51 


0.06 






0.77 


0.05 


19.10 


0.07 












1024-1239 




0.56 


0.06 


0.57 


0.07 


0.76 


0.05 


19.28 


0.07 


2.52 


0.03 








1025-1236 




0.38 


0.06 


0.35 


0.07 


0.34 


0.05 


18.48 


0.12 


1.68 


0.07 








1027-1159 




0.30 


0.06 


0.37 


0.07 


0.37 


0.05 


17.95 


0.12 


1.80 


0.07 








1031-1244 




0.70 


0.06 


0.65 


0.07 


0.70 


0.05 


19.76 


0.11 


2.68 


0.07 








1032-1229 




0.44 
0.50 


0.06 
0.06 


0.41 


0.07 


0.45 


0.05 
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19.04 
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0.08 


1.96 
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0.62 
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0.06 










19.45 
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3 


04 
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0.28 
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1.52 
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2.04 


0.05 
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0.59 
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0.67 
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19.31 


0.07 


2.56 
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2 
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0.03 
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0.06 
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0.07 
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0.21 


2 
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0.39 
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0.42 


0.07 


0.37 


0.05 


18.49 
18.00 


0.12 
0.12 


2.04 


0.05 








1138-1142 


0.45 


0.45 


0.06 


0.51 


0.07 


0.52 


0.05 


18.81 


0.08 


2.36 


0.03 


2 


72 


0.04 






0.46 


0.06 


0.54 


0.07 






18.86 


0.08 


2.44 


0.03 








1138-1225 




0.36 


0.06 


0.30 


0.07 


0.36 


0.05 


18.34 


0.12 


1.44 


0.07 








1138-1228 




0.66 


0.06 


0.48 


0.07 






19.63 


0.11 


2.24 


0.05 








1139-1154 




0.44 


0.06 


0.41 


0.07 






18.76 


0.08 


1.96 


0.07 








1139-1217 


0.48 


0.44 


0.06 


0.49 


0.07 






18.70 


0.12 


2.28 


0.05 


2 


72 


0.04 


1145-1155 




0.52 


0.06 


0.47 


0.07 


0.44 


0.05 


19.13 


0.07 


2.20 


0.05 












0.50 


0.06 


0.45 


0.07 






19.04 


0.08 


2.16 


0.05 








1147-1252 


0.58 


0.53 


0.06 


0.59 


0.07 


0.62 


0.05 


19.05 


0.08 


2.56 


0.07 


2 


92 


0.03 






0.49 


0.06 


0.69 


0.07 






18.93 


0.08 


2.72 


0.03 


3 


00 


0.07 


1149-1159 




0.32 


0.06 


0.31 


0.07 


0.23 


0.05 


18.09 


0.12 


1.48 


0.07 








1149-1246 


0.48 


0.45 


0.06 


0.47 


0.07 


0.48 


0.05 


18.85 


0.08 


2.20 


0.05 


2 


72 


0.04 






0.52 


0.06 


0.51 


0.07 






19.15 


0.07 


2.36 


0.03 


2 


76 


0.04 


1208-1151 




0.61 


0.06 










19.47 


0.11 













34 



Cluster Galaxy Evolution and the LCDCS 



Table 2 
Redshift Data: Continued 
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0.06 











55 
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07 
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0.34 


0.06 














18.20 





12 


1230+4621 


0.51 
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53 


0.07 








19.12 





07 
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0.41 


0.06 





43 


0.07 








18.60 





12 






0.49 


0.06 





43 


0.07 








18.99 





08 


1329-1256 




0.31 


0.06 





37 


0.07 





37 


0.05 


18.00 





12 






0.34 


0.06 





36 


0.07 








18.20 





12 


1333-1237 




0.99 


0.06 





77 


0.07 








20.50 





21 


1404-1216 


0.39 


0.34 


0.06 





35 


0.07 





47 


0.05 


18.28 





12 


1405-1147 




0.96 


0.06 











64 


0.05 


20.42 





21 


1406-1232 




0.96 


0.06 














20.43 





21 


1408-1209 




0.34 


0.06 











48 
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18.19 





12 
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0.06 














19.43 





07 
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0.60 


0.06 











55 


0.05 


19.43 





07 
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0.66 
0.75 


0.06 
0.06 











71 


0.05 


19.63 
19.90 






11 

14 


1412-1222 
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0.31 
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0.06 
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12 
12 
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0.58 
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0.69 
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Nelson et al. 



Table 3 
Evolution Parametrizations 



Type 


Functional 
Form 




Cocmcients 


Luminosity 


(M} + 51og/i) = Ci + C 2 z 


Ci 


= -21.70 ±0.04 






c 2 


= -0.94 ±0.09 


Color (V - I) 


(V-I) = Ci+ C 2 z + C 3 z 2 


Ci 


= -0.02 ±0.11 






c 2 


= 6.73 ±0.41 






c 3 


= -4.09 ±0.35 


Color (I - K 1 ) 


{I-K') = Ci + C 2 z + C3Z 2 + C4Z 3 


Cr 


= -0.50 ±0.76 






c 2 


= 15.1 ±4.1 






c 3 


= -23.2 ±7.2 






c 4 


= 12.6 ±4.2 



